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Abstract 

This paper proposes a Monte Carlo technique for pricing the for- 
ward yield to maturity, when the volatility of the zero-coupon bond 
is known. We make the assumption of deterministic default intensity 
(Hazard Rate Function). We make no assumption on the volatility of 
the yield. We actually calculate the initial value of the forward yield, 
we calculate the volatility of the yield, and we write the diffusion of 
the yield. 

As direct application we price options on Constant Maturity Trea- 
sury (CMT) in the Hull and White Model for the short interest rate. 
Tests results with Caps and Floors on 10 years constant maturity trea- 
sury (CMT10) are satisfactory. This work can also be used for pricing 
options on bonds or forward bonds. 
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Introduction 



The way most practitioners use to price CMT is to consider the CMT as a 
simple function of the CMS. Then the function's parameters are calibrated 
with the spreads between the forward CMS and the forward CMT. This 
spread has been increasing from beginning 2010, until mid-2011 when it 
stopped quoting. See Figure [T] below. 
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Figure 1: historical quotations of the spread between the JPY forward CMS10 
and the JPY forward CMT10. This figure shows the explosion of the spread 
from the beginning 2010 

But now there is no more liquidity for these spreads. As a consequence 
it has become very difficult to consistently price CMT and options on CMT, 
even with the inconsistent method described in the lines above. 

The difficulty for pricing CMT properly relies in the ability of modelling 
the forward yield to maturity of the related bond. The dynamics of the 
forward yield is known. BENHAMOU (2000) writes the dynamics of the 
yield. However he did not propose how to calculate or calibrate the volatility 
of the yield. In this paper we make no assumption on the volatility of the 
yield. We actually compute it as well as we compute the volatility of the 
bond and the survival probability distribution. 

In the first section of this paper we give notations, and we remind some 
definitions. In section 2 we write the dynamics of the forward yield to ma- 
turity, through the mathematical relation between the yield and the Bond 
price. In the third section we use martingale condition to derive a partial 
derivative equation (pde) of which the hazard rate function is the solution. 
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After solving the pde, we easily compute the volatility of the bond, then the 
volatility of the yield. In section 4 we make the assumption of a constant 
hazard rate function (with respect to the time between the forward start date 
and the maturity date), and we propose a Monte Carlo routine for pricing 
the expectation of the terminal yield to maturity, and the expectation of any 
payoff on the yield to maturity. Finally in section 5, in the hypothesis of a 
Hull and White model for the short rate, we give some tests results on the 
JPY 10Y maturity constant Treasury (CMT10), and options (Caplets, floor- 
lets, Caps, Floors) on the CMT10. We draw the distribution of the CMT10 
and the distribution of the volatility of the CMT10 for various expiries. We 
also compute the Black implied volatility related to the prices of Caps and 
Floors on the CMT10. As extension, we notice that this framework could be 
used for pricing options on bonds, or forward bonds, without further deve- 
lopments. 
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1 Notations and Definitions 

In this section we give notations and we recall some definitions. 

1.1 Notations 

• B tjT = Bond(t,T,T + 9) : is the value at time t, of the T-forward 
9- Years constant maturity Bond price. We will refer to T as expiry, 
and T+9 will be the maturity; 

• y t ,r '■ is the T-forward 9- Years constant maturity yield to maturity, 
related to the previous bond; 

• Q : is the value of coupon, expressed as a percentage of the notional, 
paid by the bond at time 

• c: is the coupon rate such that the value the coupon is equal to the 
coupon rate times the time step between the last payment date and 
the now payment date: q = c(Tj — 

• CMT(T, T + 9) is the value, at time T, of c such that the value of the 
bond at time T is at par: B TjT = 1; 

• k : is the number of times coupons are paid by the bond per year; 

• R : is the recovery rate of the bond issuer; 

• DF(t,T,Ti): is the value at time t of the T-forward T { — T maturity 
discount factor: DF(t, T, Tj) = 

• Pt,T,u '■ i s the value at time t of the T-forward U-T maturity Zero- 
coupon Bond: P t ,r,u = jsr^ 

• S(t, T, Tj) : is the T-forward T; — T maturity survival probability 

• r(t): is the first time after the time t, when the bond is subjected to a 
credit event; 

• X(t, T, U) : is the T-forward U-T maturity default intensity, also called 
hazard rate function. It will be properly defined in the next subsection. 
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1.2 Some Recalls 

• Relation between the bond and the yield to maturity 



k6 

+ yt,T)« ' ( 1 + yt; 



= T,7—— 

i=l I 1 



fivt, 



With 



k8 

Ci 



And 



lei; 



In particular for constant coupons: q = c V z G (1...k6), we have 

/(x) = c 1 ~ (1 + 1 x) " g + (l + x)- fl 
(l + z)» - 1 

The inverse of the function f above will be denoted g: 

g = r x 

• The Constant Maturity Treasury (CMT) 

It is the value of the coupon rate (c) such that the Bond is at par, on the 
expiry date: Bt,t — 1- It is always defined with a constant rolling ^-maturity 
Bond. Using the fact that ATj = i; from equation we get that c , which 
is here equal to the CMT, should be the solution of the following equation: 

c (l-(l + y T , r )- g ) , 

1 = — t; rr — — + (1 + Vt,t) 

K {l + yT,T) K ~ 1 

We find that 

CMT(T,T + 8) = K (^l + y T>T )^-lj (2) 
- Vt,t 

For the CMT 9 modelling, the expression of the corresponding forward 
Bond price can be rewritten as follow 
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B t ,T = ((l + ^r,r)"-l) \ ( 1 + y f \ + (l + y t , T 



1 + yt,r) K - 1 

= f(y t , T ) (3) 

With 

/ : x i > ((1 + Vt,t) " - l) + (1 + *)"' 

7 (1 + rr) « — 1 



And 



This last version of function / will be used when pricing the CMT9 or 
options on the CMT9. 

• Hazard Rate Function (or default intensity) 

The hazard rate function is defined as follow: 

= Um F l{ T<r( t ><T + AT/r {t )>T) 

y ' AT^0+ AT 

P t (T < r(t) <T + AT) 
at-^o+ AT F t (r(t) > T) 

S(t,T + AT) - S(t,T) 
atTo+ AT S(t,T) 
1 <9 2 S, 



S(t,T) 



(*,T) 



<9t> 



-(t,T) 



And given the fact that S(t,t) = 1 : no default has occured at initial 
time, we get that 

S(t, T) = e -ljMt,v)dv 
Similarly, we express the forward hazard rate as follow 

X(t,T,U) = E? T [X(T,U)] 

And 

S(t, T, U) = e -S u T Kt,T,u)du 
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• Zero-coupon Bond 



The zero-coupon bond is the value of a contract that pays 1 at matu- 
rity. Under the risk-neutral measure, the dynamics of the zero-coupon 
Bond can be written as following: 



dP t ,T 

Pt,T 



r t dt + <j P (t,T)dW t 



— <jp(t,T) is the volatility of the zero coupon bond. In this frame- 
work we will suppose it to be known; 

— Pt t = 1 

The T-Forward zero coupon Bond of maturity U is defined as follow 

Pt,u 



Pt,T,u - 



t,T 



Then we have 

dPt,T,u 



-S(t,T,U)dt + <T P (t,T,U)dW t (4) 
With 



^(t, T, U) = a P (t, Tf - a P (t, T)a P (t, U) 

And 

<r P (t,T, U) = <T P (t,U) - a P (t,T) 

2 Yield To Maturity dynamics 

Since the spot bond (B tjt ) is a tradable and replicable asset, it drifts at the 
risk-free rate under the risk neutral probability. Under this measure, the 
dynamics of the Bond can be written as follow: 



dB t7t 



= r t dt + a B (t)dW t 



B t ,t 

Where the process (W t )(t>o) is a Brownian motion under the risk neutral 
probability. 

And since B t T is the forward, it should drift at zero, under the same 
measure: the forward is a local martingale under the risk neutral probability. 
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^ = a B (t,T)dW t (5) 

Proposition 1 

The volatility of the forward yield can be expressed as a function of the volatil- 
ity of the forward bond through the following formula 

*v(t,T)= fi ^: T \ a B (t,T) (6) 

yt,Tj {yt,T) 

The dynamics of the forward yield under the risk-neutral measure is the 
following 

dyt,T ^yt,rf (yt,r) 2/. m\rj i / * r7~i\ jt t t 

= ~o^v ^o- {t,T)dt + a y (t,T)dW t (7) 

yt,T * j [yt,T) 

Proof 

From equation <Q or ([3]) above, we can write that 

Vt,T = / _1 (5 t ,r) 

= 9(B t , T ) 

Using Ito'o lemma, we thus get 

dy t ,T = g'{B t)T )dB t)T + -g" (B ttT )d<B. tT ; B. iT >t 
Using equation ([5]), we get that 

= o m \ B t,T a B[t,T)dt+ — a B (t,T)dW t 8 

yt,T 2 g[B tiT ) g[B t ,T) 

We remind that 

9(b) = r\b) 
9(b) ' 



f"(f- l (b)) 



g(b) = 

[ftf-Hb))] 

Replacing in equation (|8j), we get the two results: □ 

dy t>T _ ly t>T f"(y t ,T) ( f(y t ,r) _ Mt;T) V dt+ fJM {tjT)dWt 



yt,r 2 f'(y t ,T) \yt,Tf'(yt,T) ' / yt,Tf'(yt,r) 
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We remind that our objective is to be able to diffuse the forward yield to 
maturity. At this stage, we have written the yield's dynamics, but our goal 
is not yet reached, since we don't have the volatility of the yield, and also we 
don't have the initial value of the forward yield. 

Proposition 2 

The initial value of the forward yield to maturity is the following: 



, 1 - S(0,T,T + 9)DF(0,T,T + 9) - R*T(0,T,T + 9) 

VO T = a 1" 1 

1 ZtiS^T^DF^T,^) 



With 



r(t,T,U) = J U P t)T)U d u (l-S(t,T,u)) 



u 

Pt,T,uS(t, T, u)X(t, T, u)du 

Jt 
Proof 

Given the initial rate curve, the initial default probability distribution and the 
recovery rate of the issuer, one can express the initial value of any constant 
maturity forward Bond as following: 



k8 

5 ,t = S(0, T, T+6)DF(0, T, T+9)+J2 K?i ~ T^S^, T, T i ),DF(0, T, T i )+i?*r(0, T, T+9) 

i=i 

In particular, for the CMT9 related Bond, we have that (equation ([2])) 
c = K ^(1 + Vt,t)^ — lj; an d by definition this bond should be at par. We 
thus find the following equality (with (Tj — Tj_i) = -): 

S{0, T, T+9)DF{0, T, T+6)+(fl + yr?)* - l) J] S{0, T, T i ) J DF(0, T, T,)+ J R*r(0, T, T+9) = 1 

i=i 

And we have the result. □ 

yo t T is a first good approximation of the forward yield. It is the value without 
the correction due to the convexity (convexity adjustment). 
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3 Calculation of the volatility of the forward 
Bond 



In this section we calculate the volatility of the forward bond, when the 
volatility of the zero coupon bond is known. And also we suppose a deter- 
ministic Hazard rate function, with respect to the time between today and 
the expiry (forward start) date. The bond's volatility will be used for com- 
puting the volatility of the forward yield to maturity, through equation (|6]) 
above. 



Proposition 3 

Under the double conditions that the volatility of the zero- coupon bond is 
known, and that the Hazard Rate Function is deterministic with respect to 
the time between today and the forward start date, we have the following 
expression for the volatility of the constant maturity forward Bond: 



<TB(t,T) = -^—(e-fT +e ^T,v) d vp ap{t ^ T + 

f{yt,T) 



c 



+R 



T+e 



i=i 

e~ ^ Xit ' T ' v)dv P t ,T,u<7p{t, T, u)X(t, T, u)} du) 



(9) 



And X is the solution of the following pde: 



R 



X(t,T,u) 



u dX 
. di 



{t,T,v)dv + £(t,T,u) 



dX 
di 



(t,T,u)^P t , T , u e-fTMt,T,v) d v du 



T+e 



dX 
di 



(t, T, v)dv + f (t, T,T + 6)\ P t ,T,T+ee- ^ 



- +a \(t,T,v)di 



(10) 



+ l E (- 1 f t (*. T ' + T ' T *)) p ^ e ~ !? x{t,T,v)dv = 
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Proof 

Given the forward survival probability distribution and the forward zero 
coupon bonds, the expression of the constant maturity forward Bond is the 
following: 



_, K0 

c 



B t , T = S(t, T, T + 6)P t ^ T+e + - V S(t, T, Ti)P tTT + R * T(t, T,T + 6) 

K *■ — ' 

i=l 

Differentiating with respect to the first variable (£), we get 

_ k6 



i=l 



+R 



e ^ x(t > T ' v)dv P t ,T,uVp(t, T, u)X(t, T, u)\ du)dW t 

fl'+O / r ru ai -| ai \ 

(Rj^ (\(t,T,u) - J ^(t,T,v)dv + £(t,T,u) + ( ^(t,T,u))P t ^ u e-ttHt,T,v) d v du 



K ^ 



r 8t 



(t,T,v)dv + £(t,T,T 



*t,T,T+ee JT 



dt 



(t,T,v)dv + ^t,T,X i ) P t:T:Ti e 



■fP \(t,T,v)dv 



)dt 



with 



S(t,T, U) = e~JT 



X(t,T,v)dv 



Then the volatility of the forward bond is the Brownian motion coefficient 
divided by the value of the forward bond. 

Since the forward bond should be a local martingale, the drift term of 
the expression above should vanish. We thus get the pde, and this ends the 
proof. □ 



Looking at equation (10), we see that the hazard rate function and thus 



the survival or default probability distribution of the bond issuer depends on 
the coupon rate of the Bond. This is a limitation of this framework, and this 
observation shows that the hypothesis of a deterministic hazard rate function 
is obviously wrong. To be perfect one should use a stochastic hazard rate 
function, and this will include the volatility of the time-to-default, and the 
correlation between the time-to-default and the interest rate. 
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Now we do have the dynamics of the forward yield; we have the initial 
value of the forward yield; we have the volatility of the forward bond and 
the volatility of the forward yield through equation ([6]) . We can thus diffuse 
the forward yield to maturity and price any derivative on this underlying. 



However, it is not easy to solve equation (10), and when dealing with 
the constant maturity treasury, the coupon rate is a function of the terminal 
value of the yield, and we don't have it. 

4 Simplification of the hazard rate function and 
proposed algorithm for diffusing the yield when 
dealing with the CMT 

In this section we make the assumption that the Hazard rate function does 
not depend on the time between the expiry (forward start date) and the 
bond's maturity. The coupon rate here is the CMT9. In this context we 



propose a numerical solution for equation (10), and finally we propose a 
consistent algorithm for diffusing the forward yield, and to price, in Monte 
Carlo, any pay off on the CMT9 



We actually write that 
Then we have 



\(t,T,U) = \(t,T) 
S(t,T,U) = e-( u ~ T ^ 



Equation (10) becomes: 



n / T+ [xd-.T) 

T 



-{u-T)°^{t,T) + i{t,T,u) 



+ ((1 + Ht.t)= - l) E H T * " ^(t,T) +£(t,T,T 4 ) J P t> 
i=i ^ ' 



-(Ti—T)\(t,T) = q 



Still the pde is not easy to integrate. In the next lines we propose a 
numerical resolution. Given a discrete subdivision of the time space between 
and T, we use the following approximation for the first derivative: 

ax AfwT-Afe^ 

dt Kj J 2Atj 3 ~ 
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We input this in the pde and we get the following expression for the 
hazard rate function at each time step: 



\(t j+l ,T) 



^(tj-i,tj,T) 



With 



+ 1(1 + y T ,r, 



jf (A(/ r :/-) 



(«-T) 



e -(u-T)X( tj ,T)p t Tudu 



*(*;, r) = ^ eAfe ' T) ^,T,T +e + ((1 + Mt,t)= - l) E ^^^e-^^^P^^ 



T+6» 



[1 + {u - T)\(tj, T)} e-^- T ^^P t ^ u du 



• X(t ,T) = — |ln f S g^x)^ ) * s an estimation of the initial forward ha- 
zard rate; 

• A(t_i,T) = 0; 

• = — tj—\. 

We have everything now to propose the final algorithm for the pricing. 



Algorithm 

1. Compute the initial forward yield yo t T, as specified in the proposition 

2 above; 

At each time step: < tj < T 

2 Diffuse equation (|4]), for all maturities from T to T + #: T < u < T + 9; 

3 Compute the Hazard rate function, and the survival probability as in 



expression (10). Use yt _ 1} T f° r computing the coupon rate: y^.r 

ytj-uT', 
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4 Compute the volatility of the forward bond <TB{tj,T), using equation 
(|9]). Use y tj _ u T for computing the coupon rate: y T)T = y tj _ lt T] 

5 Compute the volatility of the yield using equation <^ : 

°« (ti ' T) = y^Mbif B(ti - uT) 

6 Diffuse the yield to maturity using equation ([7]) , and move to the next 
step; 

7 At maturity calculate de CMT using equation (|2|). 



5 Application and results in Hull and White 
Model 

Hull and White model 

In this model, one makes the assumption that the short term interest rate is 
a mean reverting, normal distributed process: 

dr t = a(f — r t )dt + adW t 

• a is the volatility parameter. It is supposed to be constant; 

• a is the mean reversion parameter. It characterizes the speed of rever- 
sion of the interest rate toward r. the bigger is a, the faster the short 
rate reverts towards f. It is supposed to be constant. 

• f is the level toward reverts the short interest rate; 

f could be deterministic or stochastic(2 factors). When f is not stochastic, 
then the volatility of a zero coupon Bond with maturity T is given by: 

a 

And the volatility of the forward zero coupon Bond is 

a P (t,T,U) = -(e- a F- i )-e- a <P- t '>) 
a 

Once we have calibrated Hull and White model parameters on market ins- 
trument prices, we can calculate the CMT using the algorithm described in 
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the previous section; and finally we can price derivatives on the CMT6. 



Tests and results 

All the tests below have been done in JPY currency. The initial hazard rate 
function (piece-wise constant) has been stripped from the spot prices of the 
Japanese Government Bonds (JGB). The tests have been realized on the 
same date: 28/03/2012. The rolling constant maturity is ten years: 9 = 10 
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Figure 2: test of convergence. 1 year expiry: 28/03/2013; R = 20%; dis- 
cretization step = 1 day; a = 10%; a = 1%. This figure shows that the 
price of the terminal yield converges very quickly: from 512 paths we have 
satisfactory price. 



ATMF Caplet price convergence 
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Figure 3: test of convergence. 1 year expiry: 28/03/2013; R = 20%; dis- 
cretization step = 1 day; a = 10%; a = 1%. This figure shows a rapid 
convergence for the ATMF (at the money forward) Caplet price: from 128 
paths. 
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sensitivity to the HW volatility parameter 
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Figure 4: sensitivity to model Parameters. Expiry = 28/03/2013; Matu- 
rity =23/06/2013; R = 20%; discretization step = 1 day; a = 10%; Mtcl 
Paths=1024- On this figure we see that the expectation of the terminal yield 
decreases with the short rate volatility parameter, whereas the ATMF caplet 
price increases. 
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Figure 5: sensitivity to model Parameters. Expiry = 28/03/2013; Matu- 
rity (Caplet) =23/06/2013; R = 20%; discretization step = 1 day; a = 1%; 
Mtcl Paths=1024- On this figure we see that the expectation of the terminal 
yield increases with the short rate volatility parameter, whereas the ATMF 
caplet price decreases. 
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sensitivity to the recovery rate of the issuer 

1 :».o 




Figure 6: sensitivity to Recovery rate. Expiry = 28/03/2013; Maturity 
(Caplet) =23/06/2013; discretization step = 1 day; a = 1%; a = 10%; Mtcl 
Paths=1024- On this figure we see that the expectation of the terminal yield 
decreases with the recovery rate, whereas the ATMF caplet price is almost 
constant with the recovery rate. 



Yield's Convexity Adjustment as f unction of 
the initial forward yield 
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Figure 7: sensitivity of the convexity adjustment to the initial value of the 
forward yield. The convexity adjustment being the difference between the 
expectation of the terminal yield, and the initial value of the forward yield. 
Expiry = 28/03/2017; discretization step = 1 day; a = 1%; a = 10%; R = 
20%; Mtcl Paths=1024- On this figure we see that the convexity adjustment 
is almost constant with respect to the initial forward yield. 
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Term structure of the convexity adjustment 
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Figure 8: discretization step = 1 day; a = 1%; a = 10%; Mtcl Paths=4096; 
R = 20%. We observe on this figure that the convexity adjustment increases 
with the maturity 



terminal distribution of the volatility of the 
constant 10 Years maturity yield 
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Figure 9: discretization step = 1 day; a = 1%; a = 10%; Mtcl Paths=8192; 
R = 20%. We observe on this figure that the volatility of the yield to maturity 
is log normal, or very close to a such distribution. 
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1 Year Expiry 


5 Years Expiry 


10 Years Expiry 


Min 


19,21% 


10, 00% 


9, 09% 


Max 


493, 61% 


495, 10% 


460, 21% 


Average 


77, 78% 


43, 70% 


35, 08% 


Standard deviation 


60, 28% 


39,48% 


31, 06% 


skewness 


309, 78% 


493, 99% 


485, 28% 


kurtosis 


1220, 45% 


3588, 85% 


3672, 24% 



Table 1: Statistical data for the distribution of the yield to maturity's volati- 
lity 



terminal distribution of the constant 10 Years 
maturity yield 
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Figure 10: discretization step = 1 day; a = 1%; a = 10%; Mtcl Paths=8192; 
R = 20%. We observe on this figure that the yield to maturity is not log 
normal. 





1 Year Expiry 


5 Years Expiry 


10 Years Expiry 


Min 


0,07% 


0, 10% 


0, 12% 


Max 


3, 76% 


8, 03% 


8, 70% 


Average 


1,20% 


2, 25% 


2, 77% 


Standard deviation 


0, 57% 


1,08% 


1,31% 


skewness 


45, 25% 


48, 22% 


40, 34% 


kurtosis 


-5, 80% 


-4, 87% 


-16,46% 



Table 2: Statistical data for the distribution of the yield to maturity 
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Figure 11: Black implied volatility surface, related to Caps and Floors prices 
on the CMT10. Caplets or Floorlets are paid with a frequency equal to 3 
Months. Discretization step = 1 day; a = 1%; a = 10%; Mtcl Paths=1024; 
R = 20%. 

Black Implied Volatility Smiles: short term 
maturities 



28/09/2012 

28/12/2012 

28/03/2013 



2,00% 



585,00% 
535,00% 
485,00% 
435,00% 
385,00% 
335,00% 
285,00% 
235,00% 
185,00% 
135,00% 
85,00% 




0,00% 



0,50% 



1,00% 
Strike 



1,50% 



Figure 12: Short term maturities Black implied volatility smiles, related to 
Caps and Floors prices on the CMT10. Caplets or Floorlets are paid with a 
frequency equal to 3 Months. Discretization step = 1 day; o = 1%; a = 10%; 
Mtcl Paths=1024; R = 20%. The curve is very important on the short term 
implied volatility. 
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Black Implied Volatility Smiles: mid-term 
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Figure 13: Mid-term maturities Black implied volatility smiles, related to 
Caps and Floors prices on the CMT10. Caplets or Floorlets are paid with a 
frequency equal to 3 Months. Discretization step = 1 day;o~ = 1%; a = 10%; 
Mtcl Paths=1024; R = 20%. 



Black Implied Volatility Smiles: long term 



70,00% 




10,00% 
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Figure 14: Long-term maturities Black implied volatility smiles, related to 
Caps and Floors prices on the CMT10. Caplets or Floorlets are paid with a 
frequency equal to 3 Months. Discretization step = 1 day;a = 1%; a = 10%; 
Mtcl Paths=1024; R = 20%. 
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Conclusion 



In this paper we have used the martingale results to derive the dynamics 
of the forward yield to maturity, and to calculate the volatility of the for- 
ward yield to maturity. This is possible when we have the volatility of the 
zero-coupon Bond. We made the assumption of a deterministic hazard rate 
function, with respect to the time between today and the expiry date. We 
have proposed a result for the initial value of the forward yield to maturity. 
These results enable us to diffuse the yield to maturity, and to price deriva- 
tives on the forward CMT or the forward bonds. 

As application we have supposed a Hull and white model for the short in- 
terest rate, we have computed the volatility of the zero-coupon bonds, then 
we have calculated the convexity adjustment on the yield to maturity, we 
have calculated the terminal yield to maturity distribution, the distribution 
of the terminal volatility of the yield to maturity, and Black implied volati- 
lities surface related to Caps and Floors on the CMT10. 

As results we observe that the convexity adjustment of the yield increases 
with the maturity. We observe that the volatility of the yield to maturity 
seems to be log normal, whereas the yield to maturity itself is not log nor- 
mal distributed. The black implied volatility generated is decreasing with 
respect to the maturity, for any fixed strike. The implied volatility smile has 
an important curve for short maturities (until 1 Year maturity), and then 
the curve decreases as maturity goes higher. 

As extension to this work, one should release the hypothesis of determi- 
nistic hazard rate function. This will includes the volatility of the Hazard 
rate function, and the correlation between the time to default and the interest 
rate. 
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